DinTarBUTION  STATRM’fOT  A 


Approved  for  public  te!(>096| 
Distribution  Unlimited 


jrjVAtNicf't* 


r- 


iMAUYTICAl. 


1'.^ 

Iv'  jt  7HOOG 


UNIT 


I 


' * ! I 


r>  n r- 

L.  . s...'  V » 

AUC  39  mt  ' i j 


V.  u 


t «_r  Lj 


mSTRIBUTION  STAT^;ot;A 

Appioved  for  public  relecse; 
Distribution  Unlimited 


f*  r<  >:  i"'  t\ ! ( J ’•  C'  Y i i' 

Ns... 

•i 

HEAVY  o(!'_IY  aR  I EL£C.'i  ''O-'l'".}  LE  >'l ' T it  oY 
. syra:i'!;f,  ► E'’ vr 


J a/9 


a'ABUf;  OF  cu;;iv;i;To 


-j- 


Sect.]  cn 


Titin 


Paqe 


1 

2 

3 

4 

5 

6 
Y 
8 

9 

10 

11 

12 

13 

14 

15 
IC 


InlrochicLi.on 

Derivation  of  Optiir.al  i'l-occ'osor 
li3G  of  Da^ ; ' or  Information 
Doppler  Error 
Ti^nsient  Tar£;et 
Fadinc  Tar^ot 

Date.  Storaco  Considerations 
Modified  Definition  of  Detection 
Conclusions 

Ferforuance  of  Optlrral  Irooossor  in 
Additive  IJoise 

hur.ioer  of  Mixi-d  Hypotheses 

Distribution  of 


Distribution  of  Sinrjle  Fine  Processor  Output 
Under  UVTA  Conditions 

PerfoiT.iance  of  Tnrcchold-Square  Lau  Ifocessor 
In  Additive  lioice 

Fcrforrance  of  Threshold-Cciuaro  Lav/  Processor 
Under  UVTA  Conditions 

Probability  of  Detection  vith  I>ita  Reduction 
Scheme 

Referoncv's 


1 

2 

9 

l4 

19 

21 

2y 

29 

30 

31 

35 

36 

36 

4i 

46 

50 


□ o 


-1- 


KU),'i'i}'j.K  n::c;  r.o::AH  ]).■■: i'cicTion 


1 . Introduction 


Idle  doGi£;^  of  detectors  for  single  I'inc  sonar  has  closely  foll.ov/ed 
rud'-ir  desicn  x^rinciples . hatched  fi iters,  optiinUi'ii  only  for  a staid onary, 
Gaussian  noise  channel,  a^’e  used  for  detection  althouxdi  ocean  noise  has 
been  shora  to  be  neith?r  Gaussian  nor  stationary.  Although  much,  has  been 
learned  about  the  sonar  enviromaent  in  recent  studies,  more  accurate  desi£';n 
and  analysis  of  detectors  inust  avalt  detailed  knov/ledse  of  the  ocean  F;ediiur. 
vhich  is  unavailable  a.t  xiresent.  The  radar  analysis  and  t'ne  desipn 
principles  derived  therefrom  are  thus  applied,  essentially  unchanged,  to 
siNrf.e  pine  sorsir. 


In  the  design  of  ciultiiile  ping  sonar  there  is  an  additional  coaaplicati on 

led  by  Marcum 
aracteristics  are 

assumed  to  remain  stationary  for  the  diiration  of  the  pulse  ti-ain.  . Tnis  is 
a reasonable  assumption  since  the  entire  pul.se  train  may  last  only  micro- 
seconds, But  in  the  20  to  30  seconds  between  sonar  pings  the  target  location 
end  other  characteristics  can  change  markidly.  The  problem,  then,  is  to 
detect  the  presence  and  estimate  the  location  of  a target  whose  path  of 
motion  is  subject  to  cert-aln  constraints  but.  otherwise  unknown  to  the 
obseivcr. 

In  this  report  a for.uulation  of  the  optlmun  multiple  ping  detector  is 
present ;d  v/hich  takes  Into  account  the  stationarlty  differences  between  tV.e 
sonar  problc;n  and  t!ie  ide.nl  indar  fo.  i ul  ation  by  .'.'nreu::. 


vhich  is  not  present  in  the  analogous  radar  problem  stud 

f p r.-  n n tho  radai 


',r  case  all  of  the  target's  ch 
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^ • I>-^rlval,it>n  oT  Optiirial  Iroaogsor 

I«t  N be  the  total  naabor  of  (discrete')  locations  under  surveillance 
by  the  sonar  system.  By  the  J pine  there  are  possible  paths  which 
the  tarcet  may  have  traveled  throu^jh  the  survcil.lancc  area.  The  data 
aval] able  frtrn  the  last  j pin^s  is  to  be  processed  in  an  optirial  fashion 
to  decide  between  the  followinc  l+ll'^  hypotheses; 

H(1):  the  target  is  pi-esent  v;ith  path  jfl 


H(N'^)  : the  target  is  present  with  path 

K(1+1I'^);  the  target  is  absent 


The  data  from  ping  k,  1 < k < J, 


can  be  surr.marizcd  in  the  vector 
ingle  ping  processor  output  which 


* 

th 

corresponds  to  the  n location  on  ping  k.  Tne  true  location  of  the  target 

"tVi 

at  the  time  of  the  k^  ping  under  hy-pothesis  H(m)  will  be  given  the  sy:::bol 

a(k,m). 

It  Is  desired  to  process  the  data  from  successive  pings  in  a way  which 
maximizes  the  probability  of  detection  for  a specified  average  mtnber  of 
false  alams.  It  is  well  known  (see,  for  example,  Reference  2)  that  the 
data  processor  which  accomplishes  this  objective  must  compute  the  set  of  a 
posteriori  probability  d'^nolty  functions  p(Zj,  ...,  Z^jH(m)),  m - 1,  l+N*^, 

and  compare  those  with  a threshold  c which  is  inversely  weighted  by  the 
probability  that  H(rn)  is  true.  That  is,  if 


v(h, 


z,  H(m)) 


{ H(m)J. 
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tlu.n  lt(m)  iji  ticcopicd  v'acrcas  If 


then  H(i/i)  i r.  rejected. 

Vfe  assirne  that  the  wavcfori.i  r(t)  returned  hy  the  target  concicts  of 
the  trariGmJ  tted  uarrovLand  r.icnal  with  envelope  V(t),  carrier  fjeepueney 
4)  Q and  rando.n  carrier  phrase  0 jdur.  additive,  stationary,  white  Ga.uGSian 
noise  v^(t)  of  spectral  density 


r(t)  = s(t)  + Wj^(t)  = V(t)  CCS  (fDj,  t - 0)  4 Vj^(t) 


The  cinale  ping  processor  consists  of  a bank  of  K matched  f ii ter  --  envelope 
detector  corribinations,  one  for  each  location  in  the  surveillance  area. 

Each  of  the  outputs  has  the  fora: 


r(t)  V(t)  cos 


Oi^t  dt) 


2 


+ ( 


r(t)  V(t)  Gin 


ro^t  dt)‘ 


where  T is  the  duration  of  the  transmitted  sirnal. 

8 

If  the  target  is  not  present  the  returned  w'avefonn  consists  of  noise 
elone  and  the  probability  density  function  of  X is  given  by  (see  Pefcience  6 
or  Section  13  of  tVie  Appendix  of  tliis  report,  with  Wg  etiuil  to  l)  : 


p(x) 


exp 
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vhere  , E is  t>ic  trtinsmitLecl.  sicnal  enor{;/  and  V/  is 

the  RF  bandv/idtli  of  the  receiver.  When  the  tar£;et  is  present  the  pri.>hahility 
density  function  of  X is  given  by: 

P(x)  .--g  exp  ( - ) I^(x), 

2 o- 


Tnat  is,  X has  a Rayleityi  distribution  with  2 degrees  of  frecdo';i  if  the 
target  is  absent  and  a Rice  distribution  with  2 deg.”ees  of  freedom  if  the 
target  is  present,  dhie  single  ping  processor  outputs  corresponding  to 
different  pings  and/or  different  locations  are  assu'r.ed  to  be  indepenrlont. 

The  op>tirrial  data  processor  must  comimte: 


p(^n#  •••#  ^2N*  •*•>  •>  •••!  H(ri)) 


'll#  •••#  -HI#  -21 

J N 


JN 


(I) 


Row 


'kn 


ct' 


^ exp  - 

2^ 


if  n / a(k,m) 


(2) 


2 4 

exp  ( - if  n = a(k,m) 

2 ex' 


Taat  is,  if  n corresponds  to  the  true  target  location  a(k,n)  on  ping  k 
vader  hypothesis  H(m),  m * then  X is  Rice  distributed;  otherwise 

it  is  Rayieigji  distributed.  Under  hypothesis  HR"  (target  Is  absent),  X has 
a Rayleigh  distribution  for  nl]  H locutions  and  all  J pings. 
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SubiitituLinc  (2)  in  (l),  ’.vc  hnw 


7 

I H(m)) 

Jr 

k-l 

K 

1T 

n--l 

Xkn 

— 2 exp 

C3^ 

kn 


n / a(k,in) 


exp  (-  I (X 


= exp  (■ 


■)  n 


N_ 

Tl  exp  - p 

n=l 


kn 


M^,e(k,m) 


m - 1,  ...,  r 


r 


and  for  n = HII^: 

p(Z^,  •••>  Zj  I H(tn)) 


f 


N 


T 

k-1  n=l 


2 

Xkn 

exp  2 

2 cT 


Those  factors  of  v{Z^f  ...,Z^  j H(m))  vhlch  are  the  same  for  each  hypothesi 

need  not  be  computed  since  they  will  have  no  effect  on  the  relative 

ina£7^1tude3  of  p(Zj^,  Zj  | H(ni))  for  ra  = 1,  l+II^.  Similarly,  the 

1 2 

constant  factor  exp  - cr-  niay  be  included  in  the  comparison  circuitry. 

2 


Tnus,  only  the  factor: 

if 

llil  ^o(^k,a(k,n)^ 


k,a(k,n 
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Eiaot  b?  co::;iiut/.d  by  the  d”l'i  proee'-.r.or  and  cor, ted  with  a thrcnliol.d  for 

m = 1,  . . . , 

Accc.,l  ir: 


Id  reject  Il(..,)  If:  ir 


If  the  teot  futiction  doer,  not  exceed  the  threnhold  for  any  of  tne 
hypotheses  Il(ri),  n = 1,  N"^,  then  the  t-Tcet  is  declared  to  bo  absent; 

that  is,  H(H  li'^)  is  accepted.  Decisions  based  upon  this  test  function 
will  have  the  desired  optin-.al  properties. 

Now  that  the  optimal  test  function  has  been  obtained,  any  rr.onotonlc 
function  of  it  r.ay  be  used  Just  as  well.  Fcr  example, 

_ J . 

n*  ^o^^k,a(k,m) ^ ~ ^o^^k,a(k,n) ^ 


The  decision  criterion  becomes: 

Accept  H(m)  If:  (j)  = ^ ^o^\,a(k,m)^  ^ ^ 


end  reject  H(m)  if:  r^j)  = ^ y^k,a(k,m)^  " ^ 


“7- 


llio  jath  which  the  tari^oL  fc/ilo'.-’i:  ir.  celectccl  by  an  iritci i i 
controllr-r  acoardirii;!  to  a jilan  v’trj'ch  ia  imhnown  to  the  sonar  c'lorntoi-, 
Tliat  is,  incivnenl.s  in  the  V*.r(j;.'L's  j'ath  arc  deteniiini stic  anl  unknov;n 
rattier  than  ranloi.  Tnc  naxinmi  aMount  tint  trii  tar^^et  location  can 
chancy  bctv.'ecn  piic^  i-’  clotemincd  by  tlic  velocity  and  acceleration 
] i.raitetj one  of  ttic  tar^jct.  in  the  absence  ol'  further  Infon.ati  on  the 
tar^jet  locatiori  on  the  next  pin^  has  a.  chance  of  beinc  in  any  of  t'ne 
positions  noi s'.iboidn';  its  location  on  t’no  last  pin^.  Tnc  probability 

I’r  ■^a(k,n)  j a(k-l,  m)J-  that  the  tarcet  v/ill  bo  in  location  a(k,m) 
on  pine  k w’nen  it  was  in  location  a(k-l,  m)  on  pin£;  k-1  will  in  £^oneral 
be  a function  of  both  k and  the  liypottiesis  H(i;i)  . 

Wo  asEUiie  t’nat  the  t-arcot  perforras  a random  walk  in  the  sonar 
surveillance  area.  That  is,  increments  in  the  tarcet's  path  from  pine 
ping  are  independent,  identically  distributed,  zero-moan,  vector  valued 
random  variables.  Tnis  may  bo  considered  a w'crst  case  situation  as 
regards  knowledge  of  targ-'t  behavior.  Ihcn: 


r 


= Pr 

|^a(l,m),  a(2,m),  . 

= P . 

r 

^a(  j,m)  j a(l,m),  a 

(2,m),  • • " / a(j"l 

. Pr 

|a(j-l,n)  |a(l,m). 

...,  a(j-2,m)| 

= Pr 

^a(j,m)  1 a(j-l,m)^ 

f Pr  1 

• 

• • Pr 

ll  I’_  /a(k,m) 
k-1  L 


a(k-l,n)j. 
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Substltutinc;  Inio  (3),  't-f’-  ' tcct  function  bjc-. 


If  the  velocity  of  the  tar£;ot  in  constrained  to  r.o:no  finite  value, 
say  V . the  tarc-et  ] oration  on  I'ing  k must  lie  vithin  a correepondit:" 
nmiLcr  II  of  the  (discrete)  positions  adjacent  to  its  location  on  pin"  k-1  . 
Tiic  test  function  f~*p;(j)  need  not  be  ccsiputcd  for  hypotheses  v.’hich  ai'e 
inconsistent  v/ith  the  velocity  and  acceleration  limitations  of  the  target . 
That  is,  = 0 for  certain  hypotheses  under  the  random  valk 

assumption  on  target  behavior.  In  viev/  of  the  arbitrarily  large  and 
negative  term  resulting  from  In  •^H(m  )}  f rin(j)  fo’'  these  hypo  theses 
cannot  possibly  exceed  the  threshold,  for  any  finite  values  of  the  data. 
Equivalently,  hypotheses  which  are  inconsistent  with  the  target  model  can 
be  disregarded  at  the  start.  Paths  for  which  the  target  lies  within  range 
of  Its  preceding  location  on  each  ping  will  be  called  "consistent"  peths. 

Expressions  are  derived  in  the  Appendix  (Section  lO)  for  the  average 
number  of  fal se  alarms  (J)  as  a function  of  the  number  of  pings  j and 

the  probability  of  detection  I’jj(j)  as  a function  of  j for  the  conditions 
Just  considered.  It  is  also  assumed  for  these  dorlvatior.s  that  the  terget 
remains  in  the  sonar  surveillance  area  for  all  j pings  and  that  all 
consistent  paths  are  equally  likely. 

Tne  terms  false  alan;i  and  detection  as  used  in  discvujsions  of  system 
performance  in  tliis  report,  merit  some  further  explanation.  A false  alar.n 
is  raid  to  occur  if  the  tost  function  corrosponJing  to  target  path 

m (hypothesis  li(ri))  exceeds  its  threshold  wh.eri  the  tatt^et  has  not,  in  fact, 
occupied  any  of  the  J locations  a(k,n),  k ~ 1,  ...,  J,  duririg  its  path 
throuji  trie  surveillance  area  or  is  absent  ent.'.rOy.  lim;;,  the  probability 
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/ 


I 


th'j.1  H(n:)  vill  be  crroncour'ly  acc-.jiU-d  ir,  rjqu't]  ^o  ilic  prubabLlily  th'it  the 
cmi  of  j rundc;:i  variablcc  vith  RHylci^h  d i otribul  ion  exceeds  a tbrcr.hold. 

A dr.' taction  is  said  to  occur  if  the  can  of  t’ne  j Rice  distributed  random 
variables  corrospondinc  to  the  trae  target  path  exceeds  its  thrcsiiold.  It 
Is  clear  tVut  there  arc  patiis  which  aereo  v;i th  the  true  tai^et  location 
on  certain  pin£C  and  not  on  others.  It  is  shown  in  the  Appendix  (Sectior^  ll) 
that  t'ne  number  of  these  r.ixod  hypotheses  which  are  consistent  will  be  srail 
in  comparison  to  the  total  number  of  consistent  hypotheses  and  they  ■■111  be 
nepilccted  in  the  derivations  of  expressions  for  syste;.-;  performance.  Only 
the  expressions  (l6),  (3.7),  (24),  (25)  and  (2u)  for  system  perforinance  arc 
affected  by  this  approximti on.  A modified  definition  of  detection  is 
discussed  in  Section  8 of  this  report  which  takes  these  nixed  hypotheses 
Into  acco’unt. 

3 • Use  of  Duppler  Tnfomration 

The  data  processor  can  be  modified  to  make  use  of  doppi.er  information 
when  it  is  available.  We  first  assume  that  t’nere  is  no  error  in  the  doppler 
measurement.  In  this  case  is  redefined  as  follcws : 

\ = ^\i»  ' ^kll^ 

where  D is  the  (exact)  valxio  of  dopplcr  associated  with  the  location 
on  pine  k.  For  a ti'.rtiet  whose  acceleration  and  turninp;  rate  are  constrained, 
the  values  of  dopplcr  on  successive  pinp;s  will  be  depe.nd-cnt.  For  simplicity, 
this  dependence  will  bo  assur-.ed  to  extend  back  only  to  the  last  pinp;.  Since 
the  noise  is  white,  the  d>  strlb-utlon  of  neasu.rod  dA-ppler  from  non-tarpet 
loca.tlcris  will  bo  uniforn  over  the  RF  ba.a. .Iwidi.'a  V.'  of  the  reC"Lver.  Tnat  is, 
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1 . /^'.n,, 

W (~w'^ 


if  n / c.(k,tn) 


■^<T-i 


r>  ^ 
^'k-l,n-' 


i r 


— e>  f \r  TO  ^ 
\ • */ 


vVicrc  rect>  (t)  ^ J 


J 

I 

M 

1 

0>  1 

[t] 

< i- 


is  V.’oolward's  reel  function  (kcferen.ee  2).  The  values  of  doppler 
corres^ionlinr'  to  different  locations  on  the  saree  ping  are  assiv.ved  to  he 
independent  of  each  other  and  of  tr.e  arripli tales  n = 1,  ...,  K; 

k = 1,  ...,  J.  Than,  the  joint  prchahility  density  fanction  of  t’ne  jli 
values  of  doppler  fro.n  II  locations  and  j pings  is  given  hy: 


p(^l]  > •••<  ^iri> 


III' 


jl^ 


jir 


_K_  N J 

p(Ppn'  ^jj  " P^^-.n 

n=l  n=l  k=l 


\-l,n^ 


fr 

k-l 


1 K 

Tf  w 

n=l 

n / a(k,m) 


P(P, 


k,a(k,n) 


\-l,  a(k-l,n)^ 


r 


1 .i_. 

j|.p  ^^'k,a(k,p) 


(v/) 


k-1,  a(k-l,-i)^ 


(5) 
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r 


I'fic  lact  rolat-ion  in  Eqv.tion  (5)  follow:;  frou  the  fact  that  all  rnear;urod 
valvicT.  of  will  lio  vjUiin  thj  oupjioi't  of  root 


Then:  p(Z;j^,  Zj  j li(n)) 


H(m)) 


~ ji ' * ■ ' jii  I ^ • • • ^ • * * » > • • • I ) 


exp  (- 


.1  If 

ll"  'll 

n-1 


Xj-n  ^kn 

---S  exp 2 

20^ 


ir  I (X,  ,,  0 

o'  k,a(k,r.)' 


1 4- 

^ k==l  ^^^k,a(k,m)  1 \-l,  a(k-l,in)^ 


(6) 


Once  a^ain,  those  factors  of  p(Z^.,  ...,  Zj  jH(m))  which  are  the  sar.e  for 
each  n can  be  cancelled  for  purposes  of  comparison.  Tne  test  function  then 
bo  canes: 


U p(l>i;,a(k,n)  "^-l,  n(k-l,m)) 


which  will  be  written: 

]r  i„(xp  p(D^  I 


for  simplicity. 
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I'no  tcr.t  fiincLjon  ir.  to  Lc  cc^n]  arcd  vltli  a thrcnhold  inversely  veijatod 
by  Uie  prior  probabilities  1'^  (l'(rO]  : 


iX  1 < 7 


{ll(r0] 


J_  r 


k=l 


/r  Ip  |a(k-],  n)J 


J /n  C 


(7) 


Tile  tost  function  now  contains  the  terns: 


/n  p(r^  Dj,.p)  P(\^a(k,ri)  1 Vl,  a(k-l,rn)^’ 


k^li  • * • f J 


which  woi[>it  n:ost  heavily  those  retui'ns  which  have  doppler  values  consistent 
vith  the  last  observed  doppler  and  with  the  path  being  tested.  The  location 
factor: 

J^n  k^ 

Is  large  if  location  a(k,m)  is  easily  accessible  from  location  a(k-i,  n) 
in  view  of  an^'  available  information  conceming  target  be'aavior  and  is  not 
dependent  on  the  data.  However,  system  perfonrance  can  be  improved  by 
updating  those  location  factors  with  the  incoming  dvoppler  infornation. 

Tint  in,  F'^  |”a(k,  n)  I a(k-l,  n)  1 is  replaced  by: 


^a(k,n)  a(k-l,  m)j 
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P’oi'  the  targets  uriicr  coriGidcratlon,  this  has  the  effect  of  rcducinc  th.c 
mtnhcr  of  consistent  paths^  thereby  reducing  the  Ei.vcra£;e  mrnber  of  false 
alari’is.  An  cxanplc  will  illustrate  the  operation  of  the  processor  witri 
uniform  location  factors  which  are  updated  with  exact  doppler  information. 

Suppose  that  the  tar£;et  is  travclinc  with  velocity  v alone  S-  straic'nt 
line  toward  the  sonar.  The  maximum  velocity  and  acceleration  of  the  tarcot 
are  assar.ed  knov/n  and  cquiO.  to  v^^^  and  a^,  respectively.  If  a(k-l,  m)  was 
the  di.stance  bctv;een  sonar  and  target  on  the  last  ping  and  T is  tlie  tine 
between  pings,  then,  assuming  that: 

T - V, 

the  distance  a(k,  m)  between  sonar  and  target  on  the  next  ping  must  satisfy 
a(k-l,  n)  - vT  - -2-  a(k,  n)  a(k-l,  m)  - vT  + a^ 

p 

That  is,  a(k,  m)  is  constrained  to  lie  within  an  interval  of  length  a^.T" 
when  the  velocity  (doppler)  of  the  target  is  known  exactly.  V/ithout  dopple 
Infornation  e(k,  m)  need  only  satisfy: 

a(k-l,  n)  - T ^ a(k,  m)  ^ a(k-l,  m)  v^^^  T 

or  a(k,  m)  is  constrained  to  lie  within  an  interval  of  length  2 T.  Tne 
ratio  of  the  mnibcr  of  consistent  paths  with  doppler  infors'.ation  to  the 
nuisbcr  H without  doppler  Iiiru2T.ation  is  thus: 


H , a T 

^ < 1 
K 2 ~ 


Replacjn*;,  II  by  In  Equation  (l6)  results  in  a lower  value  of 
Tfic  i)roli:il)ility  of  detecLion  ir.  thua  Miller  when  doppler  inforr.ation  is 
available  for  a civen  value  of  under  these  conditions.  More 
complicated  ta,rc,jt  behavior  including  toa-^^et  rotation  can  be  treated 
similarly. 


4 . Itoj^ile  ^_Erroi- 

Let  the  measured  doppler  now  consist  of  the  true  doppler  D.  plus 

icn 

the  measurement  error  w,  ; 

kn 

A 

D.  = D.  + V, 
kn  kn  kn 

Tile  jll  quantities  w^^,  k = 1,  j;  n = 1,  N are  assur.ed  to  be 

independent^  identically  distributed  random  variables  which  are  also 
Independent  of  the  true  doppler  The  data  vector  is  now  given  by; 


if  1 r i 

k=l  n=l 


A 

1 

^ rect  w‘ 


n ^ a(k,m) 


^ ^k,a(k,ra)  ^1,  a(l,m),  •••, 


n 
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^ ^ l:^a(K,ni)  J,  a(l^ri)’  a(k-l,m)^ 


(8) 


vhlch  will  be  vritton 


yjO*-iT  -j|- 


for  siiriplici  ty . Conditioning  on  vc  have: 

,A  ( A A . 

p(\|  ..., 

J •••^  \-i'  \-i^  Villi’  '*•'  ^ \-i 

w 

where  the  Intecaitioa  is  over  ’.7,  the  support  of  Then; 


(9) 


p(Djt|  Dp  ...,  ^ |\-1^ 

« P(»k  ^ '■k  1 Dk-l)  = I ^k-1^  * 

•=  P^\|  \-l)  * (10) 

Since  Wj^  Is  ln4oi»end;nt  of  both  Dj^  and  Substitutinc  (lO)  in  (9): 

A I ^ Vi> 


p(D, 


k-1 


A 

D, 


A 

Vi> 


dD, 


k-1 


> 
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I A A I A 

Approxini'itjjii'  | \ boconer.; 

p(?)^  D^,  ..., 

= / I r(l).  1)  J * D (wj  pfi).  ( n J d D 

/ L y^-1  ■ J I k-1  k-] 


vhic’n  U2)on  cu’irtJ iuLion  in  (8)  ylcldr, : 

K A 


' I I J [p(D,  . r(v^)]  I D ) d n,.j 

A-X  ‘'y  L J I 

Roplacinc  p(l\j^,  ...,  pCD^;!,  ...,  | 

canceling  corjnon  factors  yields  the  test  function  which  is  to  be  cor.iiared 


vlth  a threshold: 


^Jh)  J [pt^K  I \.l’>  * P(\>]  P(\.l  I Vi’  ^ \.l  J - 


r'(j)  = > i’nuxj 


JV— X 

■/"  /[p(\  I ”k-x)  * p'"k>]  p'\.i  I Vi>  ■*  ”k-i 

^ J^n  |^a(k,m)  a(k-l,  n)J-  ^ Jlx\  C 


n 


!Rie  ramlo.'.  varin^>3c  v’nor,e  coivUiion.'!],  lu-obability  d'-.-nnlty  function  is 
given  by  the  convolution: 

I ==  p(\  " \ I 


brsr.  varinnce  cnu'i.l  t,o: 

Var  (1^^  \-l)  ■*■  (■^k) 

and  nay  Ijo  appro;:  lira  ted  by  a uniforn  distribution  about  with  this 

variance.  The  variance  of  a random  variable  Y with  uniform  pi'obability 
density  function  p(y)  is  given  by: 

r2 

Var  (y)  = -yg' 

where  L is  the  support  of  p(y).  Trie  support  of  the  unifora  probability 

A I 

density  function  p(Dy  thus  given  by: 


1^12  [var  (Dj,  I + Var  (Wj^)]|  " 


CO  that; 


\-l^  ^ Y7  rect  


If  the  measurement  error  distribution  pC'^j^)  aesunod  to  be  uniform 
then  it  follows  tiiat: 


I ^-1^  " i-ect  — *^-"1 


whciai  L is  the  support  of  tiic  probabJllty  density  function  p(w,). 


08- 


I'iicri: 


j [p"’k  I r(  Vi  I \-i)  “ \-i 


'/  <r. 


5^  ■■  / I - \“i  N 

rect  -r:*-  ) (l~  ^ ^ Vl 


A li  - 

+ ~^- 
K 2 


.”i:i:A.-i  A 

jfroct  \ ^\.r  " 1,“^^ 


'D  - — 
k 2 


■where  /\  is  equal  "to  the  len£;th  of  the  intersection  of  the  inter.'als; 


i '■d  . "'d 


°k  - T ’ “k  5'  “■*  ”k-i  - 5-'  \-i  •*  2 


60  that: 


,A  A ^'d  ^ V X A A 

nln  (Dj^  - + ”2  - , L^)  if  _> 


A A ^ Kf 


A A 


/V  a " T \ 

“^"(^kO-\"“‘i“'^  Vk^j 


.Lj  + L 

min  (—5 K _ 

2 


A A 

^■\.i  • K') 


Substitutinc  in  (ll)  and  dropping  the  co:i’j:ion  term  jin  Lj  L^,  the  test 


function  becomes ; 


] 
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C) 


A. 


r„(j)  ■=  } w * 

k=3. 


/'d  Iv 
min  ( *2 


\ - ^- 


k-1 


, K)| 


£n  Tj.  j n(k-l,  m)J. 


whei'c,  from  Equation  (iS): 


(’ } 


4 [var  (Dj^  I D^_j)  + Var  (v^)]  | ^ 

Ihe  second  tcm  of  Pjjj(j)  is  undefined  for  values  of  |%  - %-lj 

greater  than  -i(L,  + L.,)  • I'his  corresponds  to  the  fact  that  paths  for 

d " 

which  the  mcaouicd  doppler  from  ping  to  ping  varies  hy  an  amount  greater 
than  are  inconsistent  with  assumed  target  limitations  and  doppler 

measurement  accuracy. 


5 • Transient  Tar get 

Up  to  this  point  it  has  been  assumed  that  the  target  remains  in  the 
sonar  surveillance  area  for  the  duration  of  the  pulse  train.  Tne  data 
processor  described  by  Equations  (4)  or  (T)  is  eptimua  for  this  case;  it 
will  produce  the  greatest  probability  of  detection  for  J pings  and  a 
specified  average  number  of  false  alaras  (J). 

In  a realistic  situition  j wil]  not  be  fixed  beforehand.  Tlie  sonar  will 
rather  transmit  pings  at  a constant  rate  for  an  indofinite  period.  Ti'.Is 
moans  that  the  threshold  C «=  C(j)  must  increase  with  tiic  nite.ber  of  ping.s  J 
In  onler  to  maintain  the  npoclflcd  average  m her  of  false  nlanis.  A 
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r 


trunr-ifiit  Var^jC-L  vhic'a  ciiLers  the  sm-vcillanco  area  after  £i  larco  nun'oei' 
of  pin^ri  has  teen  trans'p.i.ticd  vJlh  t’lius  have  little  chance  of  heinc 
detected.  In  orior  to  optinvilly  process  fae  echos  rotm-ned  fror'’.  a transient 
feii'cet  it  is  necessary  to  r>cstrict  the  attention  of  the  dat-a  processor  to 
the  last  several,  cay  H,  pinj^s.  Ascuiiin^  that  a probability  rr^ass  function 
P(  ^ ) ir.  knovai  for  the  length  of  tine  ^ that  the  tarp;et  cpcnls  in  the 
surveillance  a"ca,  M c?cn  be  chosen  to  naxirai'/.s  the  probability  of  detection 
of  the  data  ])rocesscr  for  a "vindov.’"  of  length  li  pin^s. 


Conditioninc  on  vo  have; 

OO 


p(M)  = ^ p(M  [ n p(  n 

i--o 


(13) 


vhero  1’(M  J ^ ) is  the  probability  of  detection  "iven  that  the  tar{iet  rerain: 
In  the  surveillance  area  for  jf  pings.  IJow: 


P(M|f  ) = 


Substituting  in  (13)  we  have; 


P(pO  = Pj.  |p(tO  >c(M)j  pj.  > m| 


M 

,^Pr(r(n  > c(M)^  pc;.) 

Ibo  threshold  C(m)  is  set  to  establish  a tolerable  false  alarr.i  rate  and  is 
Iriipllcltly  defined  througli  a rather  ccinplicated  expression  (16).  Tne  optinal 
window  si7,c  M thus  cannot  be  determined  analytically  but  nust  bo  obtained 
by  plotting  the  above  erq-rossJon  for  ?(.*')  from  curves  of  C(j)  vs.  J and 
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^ ^ choo'dii”  t.h-‘  M f r v/hlch  }{'■'•)  in  r;;jLxiir.i •,',«. d. 

IV.o  cMcf  dirficu.H-y  in  thin  niijirciach  is  the  unavaiJabil  ity  of  p(  J ) . It 
Ehould  be  recallct],  too^  that  data  r.tojT.^e  requi rcr.ionts  increase  exxorientially 
vlth  M and  may  dictute  that  a r.uboptir'al  value  of  I'  be  used. 

6 . nji  i!/  c t 

Tile  additive  Ce-uscian  noire  c’nannel  v/hich  has  served  as  a nodel  so  far 
In  thiis  repoi  L docs  not  adequately  repi’onent  a situation  in  v.'hicri  the  tary,-.t 
aspect  r.-.ay  vary  fro.n  pin^^  to  pinq;.  Tne  retmaied  si£;;:al  level  from  a grossly 
non-spherical  target  vill  build  up  and  decrease  as  the  ta.rget  pi-csents  now 
a broadside,  now  a head-on  aspect  to  the  sonar.  Tne  pcrfonnance  of  the 
sonar  systen  will  be  noticeably  poorer  under  these  circumstances.  A simple 
moilel  for  a target  with  uniformly  varying  aspect  is  a (two-dlrrenslonal) 
rectangle  which  is  orthogonal  to  the  horizontal  plane  and  v;hose  angle  with 
the  line  bet'ween  sonar  and  target  is  a random  variable  with  probability 
density  function  p(^ ) (see  Figure  l). 


Figure  1 


/ 


/ \ 

Tne  returned  v'ivcfor.n  i\tj  unlcr  urij.fonJy  varyir^  tar^ot  acyoct 
(U\'TA)  condi. tioiin  ir.  th,n  yiven  ty: 

'r(t)  = (s.'n  "{'(t))  s(t)  1 v^(t)  v^(t)  e(t)  + v^(  t) 

The  cTi  tmel  i.ic-lol  in  necn  to  now  Jnclulr;  the  multiyli cativc  noise 
as  VC']]  as  the  additive  noise  Wj^(t). 

Trie  randa . varia'ble  is  aosu'.  ed  to  be  chosen  from  a uniform 
distribution  over  trie  intei’vvl  [^0,  7rJ  . Values  of  are  independent 
from  pinj3  to  yin^  an  I reioain  constant  for  the  duration  of  the  sonar  pulse. 

Trie  i-andorn  jirocccs  (t)  = ( j)  is  thus  continuous  valued  with  discrete 

parameter  J.  Tne  probability  density  function  of  - sin  'j''  is  shoirn  in 
Section  12  of  trie  Appendix  to  be: 

P(^2)  = (1  - ^ 0 ^ 1. 

Tne  resulting;  sincle  pin^  processor  output  X vitVi  tarcet  absent  has  a Eayicio'i 
distribution  as  before,  but  when  the  tartlet  is  present,  X now  has  a 
probability  density  function  given  by  (see  Apipcndix,  Section  13): 


P(^)  = ^ ^ yi^(n,n)  X 


m=0  n-0 


2n+l 

exp  ^ 


where ; 

^ (n,n) 


(-1)"'  (1;  ?;  :rdn) 

m ! (n  1)^  (^^n)l 


(IM 


a nc] : 


J f tH  n : 0 


(];  ?:  r :n)  = J 

' n I n 

1 ^ P(r-3)  if  lu+n  ^ 0 

I'-l 


•f  i ^ ^ ,1  . 1 V.  ? . P ^ I 

CO-::>Ufy:  : 


^ r\-rN+'T*r'<',  1 -w  ^ •*  r*  f%  f 


P(^, 


M(r)) 


lv-1  IV  ]. 


r(x^n 


. J 

ir 

k=l 


K 

'll' 

n=l 

< 

n / a(l;,in) 


exp  - 


^ o 

K 

kn 


2 cr  * 


p(^  / \) 

h,a(k,rO^ 


Tl 

k=l 


N 

11' 

n=l 


~ ~ 2 
\n.  CXI.  - 

2-^ 


n / a(k,n) 


■\  2n+l  X,  \ 

IP---0  n-0 


1 

k=-l 


ir 

n=l 


~ ^2 

o.,  - i", 

2 c/^ 


Ci:0  n^O 


.5^ 


j5{r.\,n)  X 


2n 

k,a(k,n) 


r 


Canceling  factor.'  co:::.on  to  .all  hypotkcr-os  v:  cbt-'ii:i: 


t y > -V(-) 


Cl  0 n-0 


2ti 


k,a(k,ni) 


f 


vliidi  Vj'l.].  bo  vriui.on: 


C‘0 


Il"  1.  1.  /5(n,n 

^—1  El  b n:'0 


. ^ 2n 

> \ 


for  Eii  plicity, 

^iVinr»  -iV.rt  1 rim  >' 7 +>  >-,  rv  P m i--)  n f -7  ■?  >r  i tT  •{  e < “ r>-  : ^ r.  ♦ « 


funcflor: 


f\j)  - ^ ^..  />(':», lO 

. *”'r*  — _X  X ' 


kb],  ri'O  n-0 


CO  c<* 


I 1 1 

k=l  Ei-0  n=0 


n!(njf(u,n)!  'k 


f(\) 


• . 1 ? P 

vnich  IS  seen  to  vary  with  = ^£l . 

2:JoW 

/SJ 

Tiis  vciQitinp;  funotion  f(^)  for  U\"?A  conditions  inclules  a constant  term 
which  Is  obt-ained  by  cettinj  = 0! 

f(o) . /„  2..  ■/  ° 

ir,..0  (r>l)  2'^'^ 


This  constant  can  bo  Inclulcd  in  the  comparison  circuitry;  tliat  is; 


f’cj) . f(x2  > 


Is  equivalent  to: 


[*f  ('4)  - l’(0)J  J C - J f(G)  ^ C'  (J) 


o 

Curvcr.  C'P  j'(x)  - f(o)  vr. . X r.i  c Jn  Kirujr-  (p)  vj  t,n  c 

I'CLicCi.-tcf.  ir.c  lo,.--  • r\  f;io!i  of  ruch  of  iu-.-’:,'.-  for  wbic'f. 

f(x)  - f(o)  h'i'.  • "i'oly  r ly  V-*  aj'prov  Ir./ilc  1 Vrif.- 

horj 7.0:1  L.-.!  ax'S  ir.  chovn  Jr>  (3)  To*-  ^ s=  1).  Trial  Ir.^  t)i"- 


oplinal  I'l-oc-.  r.r.or  Tor  UPi'A  c'aj;  Lioao  lias  I'lc-  cfrccl 

(asD j C'li rj  .a  v Pu..-  of  7 ro  to)  vaa!--  n lu’-aa  f;o-.  l'.<- 

>j£piri'  (?)  iMlicc'.lo:;  t tl  the  Iriroohold  c r.’!  'Ul'T  incr 

p 

r.cVlitIvo  to  noir-'  ratio  . 


Ox  tfircGhoidi  ri(_; 
fadin’  t:irt;ot. 
■••ar;a  vilh  incx’Cj 


Ficure  ? 
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Ono  of  the  dc tec  Lion  nclieniec  which  liar;  been  proposed  for  use  with  a 
fading  tar^;jt  is  th"'  "r-out-of-j"  processor  (see  hefci'cnce  3)*  7n  this 
Echer.e  the  tarcet  is  declared  to  he  }n-csout  with  p.ith  n (hy})othesis  H(n)  is 
accepted)  if  r of  the  single  pine  procenr.or  outputs  correspond! nc  to  the 
J locations  of  path  ni  exceed  a threshold  c.  ITiis  operation  can  ho  written 
fon'ally  ns  follcn/s: 


Accc]'t  ll(r.)  if: 


I 

k=l 


u(Xj,  - c)  > r 


.L. 


and  reject  H(m)  if:  ^ - c)  < 

k-1 


vhere 


u(x)  = J 


1 If  X ^ 0 

0 if  X < 0 


Trie  wei{;>itin2  function: 


■> 


correspondin-3  to  this  schoraC  is  praphed  with  c = 2.0  in  Fi£pare  (4)  toi;cther 
with  the  optlrnl  weichtin"  fiuiction  derived  above  for  If/TA  conditions  with 


a larce  valiuj  of  the  additive  3ip;iial-to-noir.c  ratio  ( 


10). 


r / <■.  '■- 


Figure  4 

Both  processor",  arc  seen  to  give  very  little  v/eignt  to  'weak  i-e turns 
< c)  . Hovcver,  the  optimal  processor  gives  stronger  retuivis  a proportion- 
ately greater  vei£^it  wiieroas  the  r-out-of- j processor  velglats  strong 

returns  unifomly.  llie  two  processors  differ  nore  radically  for  lover 

2 

values  of  as  can  be  deduced  from  Figaire  (2).  For  j = 1 the  r-out-of-j 

processor  and  the  optimal  processor  arc  identical  since  both  processors 
merely  compare  the  single  ping  processor  output  X with  a thieshold.  It 
has  been  sho’wn  in  pi-evious  work  (Reference  8)  that  the  r-out-of-j  processor 
performs  best  when  j = 1,  From  Figui-o  7 it  is  seen  that  the  optimal  processor 
perrormance  is  enhanced  as  j is  increased. 

Expressions  for  the  perron.-ance  of  a sq,iare  law  processor  with 
threshold''d  dita  are  d-''rivel  in  Sections  l^r  and  15  of  the  Appendix  for  the 
caiditlons  of  additive  noise  and  additive  plus  multiplicative  noise  (IT.TA 
conditions) . 


> 


7.  r.Lorr_Q2. 

1-? 

Tnjrc  arc  II  extcnr.ioar.  of  the  IiH'  coar;ict'jr*t  paf-ir  on  pjii[i  j-1 
vhicli  rer.ul  L in  the  IIH'^  ^ consiotont  jetho  for  pi  nr;  J.  Tno  ini'^  ^ val  uer 

pi  ^ 

of  I (j)  for  pii'C  j can  bo  obtained  rccurrivcly  from  t>;c-  Il.-r’  vabucc  of 

r*  (j-1^  on  uinr  j-1  bv  adilinr  the  nnnrOT>r  in-t.e  term  to  r f j-1  ) . Ho’>?evci-. 

1-]  th 

Mi  ■ qunntitiec  must  still  be  stored  on  the  j pin^  and  so  t:io  storarjo 

roqu'. iCTicn'ts  increase  exponentially  wl  Lh  the  imebor  of  pinc^. 

Since  it  has  been  assu:  od  that  only  one  tar^;et  l.ics  within  the 
surveillance  area,  necessarily  all  or  all  but  one  of  f/ie  sincle  pine 
processor  outjaits  correspondinc  to  the  II  consistent  path  extensions  result 
from  noise  alone.  It  my  be  desired  to  evaluate  the  test  function  P(j) 
for  only  of  t'ne  H consistent  path  extensions  and  dunp  the  data  stored 
for  the  rciiiaining  consistent  path  extensions.  Tnis  v;ill  decrease 

the  niiaber  of  stored  quantities  after  j pines  from  Ifd'^  ^ to  Wyx-  ^ ; that 
Is,  by  a factor  of  yU-  path  extensions  are  chosen  accordlne 

to  a inaximun  likelihood  ranking  procedure.  This  naans  tiiat  the  path 
extensions  which  have  largest  test  function  /^(j)  are  retained.  Expressions 
for  tb.e  resul.ting  system  parforrance  as  a function  of  II  and  are  given  in 
Section  l6  of  the  Appendix  for  the  special  case  of  uni  for::!  clopplnr  and 
location  factors  and  for  both  the  additive  and  UVTA  channel  no-dels. 


8.  Modified  P^finltion  of  Detection 

A "detection"  has  been  defined  to  be  th-e  event  that  tie  integrated 
returns  from  the  true  target  path  exceed  a thrc;saold.  Other  definitions 
of  a d'.'tection  arc  reasonable.  For  example,  if  the  threshold  is  crossed 
by  the  returns  fro-:;i  a consistent  path  which  agrees  vi..h  tl  e true  target 
j)ath  at  J-k  points  and  differs  fro-.i  the  true  tc-rget  j ath  at  k points  at 
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vliicii  it  ir.  vdthin  r.  r.pocjTic-d  dir-tuiicc  of  tuc  triiv  tur^j'-t  Incatio.-i,  then 
Hue  r..'\y  alr.o  bv  called  a detection,  Tiiat  ie,  th:;  b.y pother, I r,  (joth)  v/hjeh 
bar.  1)0011  ficcci)tr.'d  differs  only  clifyi'ly  from  the  tnu  hy^  other  ir  (true  tercet 
path).  Kore  c.Tiei'allyj  a detection  i.uy  bo  dofir/'d  to  he  the  event  of  a threri 
crorrinc  by  the  tort  function  correrpondin^  to  any  clo.nent  of  a ruhcct  -Cl 
of  the  oricinnl  conristent  hypot’norcr.  In  t’no  example  above,  ic 

defiiicd  as  follo-.;r:  ilypot'/'Eir  ir  contained  in-vT  if  eac'i  point  of  }.'(rr,) 

Is  i.’ithin  G of  the  true  t?.r{^'et  location,  Tno  quantity  G.  mart  bo  mail 
cnou^i  to  insure  that  all  elcMentr  of  .0^  are  consistent.  /)-  has  elc;..onts 
under  this  definition  v.i;ci'c  h is  the  menber  of  (discrete)  locations  v/ithin 
radius  G of  the  true  t;ir"et  location. 


9.  union 

Tiie  for;.i  of  the  theoretically  optirnua  iita  processor  for  multiple  pine 
sonar  has  been  derived  under  certain  stcitod  aGGu..'.ptions  about  the  channel 
end  tfircet.  Suboptiiral  schenes  have  been  disciussed  and  trie  resulting 
perfoiTiance  has  been  analyzed  and  presented  in  the  Appendix  of  this  report. 
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AJ'PHlIOjX 


i 
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10.  PcrrornariCtj  of  Qpt;i;xil  iTo^essor  in  Adlitlvo 

Since  all  pc.f’ns  arc  equally  likely,  the  optimal  procennor  rnurt 
only  ccr.pute : 

J.. 


P(J)  = ^ /n  Ijx) 


k=3. 

for  all  conoiclent  hy])othecec.  For  X <<  1 the  follovins  approxiratlon 
is  valid: 

/n  Iy(x)  i-  y? 

Tnat  is,  the  optiral  processor  C''n  he  replaced  for  S'rall  X by  a square 
law  detector  (the  constcmt  factor  cot.  be  d^e^'ped  since-  it  appears  in 
tViC  i/.-st  function  of  each  hiypotheci s) . Vnis  ajiproxF'iition  is  not  valid 
for  laro’  values  of  X but  '..’ill  be  used  to  sir.pl ify  the  analysis.  It  r.ay 
be  ar£pj-'d  thiat  t'le  do£;:a'latioii  incurn  l by  thi.s  approxir.ation  is  less 
irportant  when  X is  larL’'  and  detection  is  alroa.l;,'  assured.  Tlius: 

J .1. 


P(J)  - ^ i„(x)  i Y y-‘ 

xT  k-'i 


If  noise  alo:e  is  present  X has  a F.ayleit;;'.  distribution  and  the  character- 
istic function  of  is  (fiefcrf-zice  If): 

S (u)  =. ^ "o' 

2 

Since  X Is  in.iepanlont  fro.i  pii’C  to  I'inj,  X is  also  i ixl^j er.l-'nt  fr‘:i  pi/qi 
to  pine  '--‘f  c'naractoristic  fuiction  P(j)  Is  ^iv-.-n  by; 


Tnjr  cun  lio  inverted  by  uritif;;  I’uir  cif  RoferGneo  (5)  to  obtuir.  the 
prebub i ] 1 ty  dear.ity  function  p^(')  of  P(j)  vhen  noir.o  alone  is  pro  cent; 


(T-’i)T 


If'  V 

0 exil  - p , ] > 0 


'ibe  probability  t’.nt  f (j)  execedo  the  ibronhold  0 for  a parti c.ul 


ular  riOifio 


path  i c : 


(y)  dx 


(2cx'V*^(j-*l)7 


vJ-l 

Y exp  - p dX 

2c3-'  ■ 


1 - I ( 


2c^'  VT  ' 


, J-1) 


viicrc : 


I , n-1) 


n-1  .j 

y e 


(n-1)' 


1e  I'earcon'c  Inco::;i.ilete  Garvaa  Function  Ratio. 

There  are  conniutent  jc,thr3  after  J pin^s  h.ave  been  tranonitted. 

Since  the  test  functionc  n ^ 1,  ...,  ini*^  coi  roopenJ  in^  to  t'n.osc 

patlis  are  approxliiiatoly  independent  for  Inrce  H,  the  p’o'nubllity  that  n of 


theoo  result  in  faiue  alun.'.r.  beco:n 
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n=l 


y-\ 


1 . 1 ( — q. — , j.i) 


I ( — -q.---,  ■-:)! 
2c-^-^  J ' 


Foi-  tho  truo  tcire-'t  r^th  X has  a Rice  distributicn  vith  n.cau  aad 

2 

variance  boUi  cnual  to  (Rcfei'once  6)  and  the  cl;aractcristi c function  of 

p 

X Is  (luferonco  ^)  ; 


1 

2 exp 

1 - ?iucr' 


froia  which  vj  obtairi  th.e  c'r.ar.nctei'i.stic  fnr’ctLon  of 


ru: 


exp 


1 ?iu  c- 


5 

1 - 2iu. 


I'liis  can  be  J nvor Led  by  ir.irip;  i'air  upO.O  oT  L’ofei’crice  p to  obtain  ti^e 
proba'ojiity  density  function  of  T'C  j)  for  the  true-  tnrQ-t  p'-tb: 

p 

exp  - .1  <r-r  /p  j,  y /1“1  y" 

v,m-  ( ---p)  (YT7M, 


'tVi 

vhoTC  Tn(')  is  tho  r.x'-lificd  B'jcgcI  functiou  of  j order,  Tt 


pen : 


C**' 

Vj)  -y  Pi(y)dr 


Jfii 


exn  (-  ,i  ca ^ \ ^ ( y \ t r ' -1\  v 

. -y-  (.  , ( . _ y ) (^rj  j-')  a ) 


Ir.  2n  J 

(2cr>  " ) 


2 c- 


, P , -^C'’ 


(IT) 


vhere; 


k 

I 


«j  <"< 


th 


r 


2 

t (p)  exp  {-  2 ) (^-)  <1  t 


Is  Ivarci'".'!’  J Q functuai . llic  p.'nr.bnbf  1 Ity  of  ib.'t.'ct  ■' on  ’^(  ,1)  is  plott..’ 
vs.  the  r.vere^,,-  nu.ber  of  fair..-  ala?...;,  ‘ypi,(.-7  ’■''it'-  ,j  as  a pirn 
Figiirv  5. 


■:V  en 


~3> 


11.  Hvfic^  bcr'  ". 

A nixed  j.'iiii  (liyjioVnocir:)  ir-  cne  v'.iicli  n^rccr.  with  tiic-  true  terijot 
path  on  k piuco  and  difi’orn  from  fne  true  t'-vri^et  path  ori  J-k  plnrs  w'ncre 
1 h j-1.  There  n'>'c  ^ conniotcat  j.^itha  in  al].  and  (U-] ) 

of  tiierc-  have  no  j.ointn  in  co:-non  vith  t’nc  one  ti-ue  tar(;et  path.  Tnc 
rcneiinini;  ir.i*^  - (II-l)  - 1 patln;  have  one  or  nore  pointc  in 

cornea  vith  the  true  t-arjjet  path. 

I<*t  li  V'C  the  nuiher  of  ci.ixcd  patiis.  Tnca: 

HI 


I(il 


J-1 


KH 


J-1 


‘ 1 - ir  <-k ) 


For  lar(;e  H: 


,H-1  j-1 


1 

1 - H 


60  that: 


N:n 


Ml 


J-1 


lL-1 

1 - N '>3.  - ~ jj 


since  II  is  larce.  In  this  case  j is  tlie  nu:ibcr  of  pin^s  beinp,  processed  and 
is  the  6a.no  as  the  lonipth  of  the  window  M of  Section  5-  Data  storac- 
re  quire  men  to  will  Unit  j to  a valu;  of  10  or  less  and  II  will  typically  be 
equal  to  100  or  nore . Thus : 

N 

_E 

Ml' 


-Vt 


r 


-36- 


12.  DirLriliu'J.O!!  of  \r^ 

Iict  G(')  110  th*j  dintribuLiori  function  of  nnd  ?’(•")  the  cli  ntrihuti  on 
fvinctlon  of  . I'no  rarido.-n  variable  ''j'  is  unifaj-;.;ly  distributed  on  [^G,'/>j 
and  V r;  sin  'j'  . 


G(v)  = 


'2-^ 

= P^  |*sin  f < 

-]  '] 

;in  V 

- + }’  ii'  > 'ii' 

J 

r)  + 1 - 

F ( '70  - sin  ^ w) , 

Ibc  probability  density  function  of  w is  then  oy: 


= d v7 


2 2 -} 

\ (1  - Wg  ) , 0 < Vg  < 1. 


13.  Distribution  of  Single  Ping  Processor  Output  und:: r Conditions 

Tnc  returned  sin"le  pin^  'ravefomi  is  civen  by: 


h 

:b  n 


7 (t)  = Wg  S(t,  9)  + Wj^(t)  • (18; 

Let  R bo  a vector  whoso  cor.ipononts  ere  the  tir.o  sa.nple  values  of  r (t;. 
Sinllarly,  let  0(0)  and  V,'^  be  the  tir.e  sonple  vectors  of  s(9,t)  nn.l  v^j_(t). 
Ijct  K be  the  covariance  n’.ntrix  of  the  tice  sateplcs  o^  f.;c  nidi -ive  noise 
Wj^(t).  Tnen  (].8)  can  be  v/ritte:.  in  t'ue  foi-p.: 

R G’(^)  ■’  ''2^ 


s(t,0)  ...  V(t)  CO.;  - O) 


COG  6-> 


G -i  I^V(t)  Gin 


sin  0 


= V (t)  COG  0 H V^(t)  Din  0 


LcL  V^.  anl  V^,  doiiotc  t’nc  vectors  of  tine  sarapied  values  of  Vj,(t)  and  V\.(t), 
rcGfectivoly . llic  single  pin"  piMconsor  consiDts  of  a rra,tc'"e.d  filter 
follov.’od  by  an  envelope  detector.  Vfnen  the  input  to  ti;o  single  pin" 

/s./  ^ 

processor  is  R,  the  outijut  X i.s  civen  by: 


V-)^  + (l/  K"^ 


where  ( ) denotes  trtinsposc. 

Tiie  quantities  R K”  and  R arc  linear  conbinations  of 

Gaussian  ra.idon  variables  and  are  thci'oforc  Gau.ssi.a,n.  From  Eqvsition  (I9) 
we  see  that  they  are  ^0  decrees  out  of  pha.se  and  arc  thus  uncorrelated. 
Since  tiicy  are  Gaussian  randor.i  variables,  they  are  fn.erefore  also 
independent . 

When  the  tarcet  echo  is  present  R = W2S(0)  -i  and  in  this  case,  for  a 
Civen  value  of  v^,  we  have; 

E (rV^  Vc)  - E(w2  S(o)  + Vp/  K"^  V^, 


Vg  s'(0)  K*’- 


BTld 


K(uV’  V,,)  = r/(o)  Vg 


Wicn  Q ~ 0 doc^rcca,  S(q)  - V an.i: 


S'(0)  K"’  V^. 


= Vg 

w;icro  cr"^  - V^,  K V^,  " ’“^s  ^ ^ Vg . Since  the  noicc  is  white  with 
Si'cctral  density  I.'^  ve  have: 


K 


-1 


Ho  W 


vherc  W is  the  RP  bandwidth  of  the  receiver. 


Tncn: 


2 <\ 
cr'  ^ 


<Vs 


N_  V/  II  W 


2 r-o  5^ 


1 

? 

I 


vhcre  E is  the  enercy  of  the  transmitted  sicna]  . 

Also: 

Vg  s'(0)  ir^  Vg  = Vg  v'  k“’  Vg  = 0 

since  is  orthoconal  to  (see  Equation  (19)  ).  V/lier.  0 = 90  de(;rees, 
S(o)  =:  Vg  and: 


-1 


*2  S (0)  K-"  . Vg  K V,  V Vg^ 


,2 


and; 


Vg  s'(0)  k'^  Vo  = Vg  Vg  K'l  Vo  = C 


-3y- 


Tnu;'.,  the  r.can'5  are  in-oporli o.iaT  to  and  90  decroea  out  of  phase, 

llie  variance  o"  R K“^  V^,  with  tari;;et  present  Is: 

Var  (rV^  V^)  = 


r-.  L’  (k'k"^  V - E R^K"-"  V )"  (R V’  V„  - E rV^'  V ) 


E (v;^  k"^  v/  v4  k"’  v^)  = 


Var  (R  K ■*  V ) 

E 


Tjic  square  root  of  the  Du.i  of  t'na  squares  of  tvo  Gaussian  rando:.i 
variables  vith  equal  variance  and  means  eqiual  to  a sin  ^ and  a cos  0, 
respectively,  has  a Rice  distribution  (Reference  7)  ’'Tith  probability 
density  function  equ-al  to: 


2 4 

V + a 


” w/w  ' ■ tvar(j.))  a (fer  'ra) 

/%✓ 

the  probability  density  function  of  X is  thus  given  by: 


(X  1 wg)  = 


~ y2  . / , 2%2 

^-5  exp  ( - 1-h  ( -2  e. 

cv^  2 C3" 


) 


~ 2 4 

exp  (- V-—  ) I w,  X) 


2 C7- 


n 


->I0- 


and  ro: 


j I '’'2^ 


d w 


'tt 


X . x2  + v/ 

~ 2 (~  2~ 
cr'  2 


_.l. 

) Iq  (Wg  X)  (1  - v^  ')  d 


• o 


Or>  o:> 


= Z I 


m=0  nK) 


m 


I (nj)^  (mui)<  2^ 


,2n-i3n  ■ "(X  ■ exp 


2 c7" 


vhcre ; 


r 


If  ll^^n  = 0 


(X;  2;  »►„)  ■’  "(  ff  ^ ^ If  ^ 0 


r=l 


Iicttinc: 


(-1)”  (I,  2;  «„) 

/(“  ^n)  = nj  (njj^  (^l^n)|  + 3n 


this  bo come c 


p(x) 


= I >: 


m-O  n-0 


2nH 

^ (n,n)  X 


exp 


. X 


2 


(20) 
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l4.  rcrTorinanco  of  'I'iu-i':;'nold-r>']_i.r.:.r:‘  I''.'..’  I'l'oc'.-r.r.or  j ii  Additive  Noi  :;c- 
I'ni:-  Eiiic3c  I'itic  procccuor  output  X is  thrcrdioldod  at  c;  that  is, 
only  those  output;;,  which  exceed  c are  lassed  on  for  nultiple  pin"  procoss- 
Inc.  I'ne  prolahJ]  ity  density  function  of  the  threcholded  random  Vc'.riabJc 

A 

X Is  pivcri  by: 


P (X)  - I 


^ (X)  + P (X)  u (X  - c) 


vl'ero  p(‘)  is  the  probability  density  function  of  X before  thresholdir; 
^(•)  is  the  Dirac  delta  function  and: 


u (y)  = 


Lot  C'(*)  be  the  distribution  function  of  X'".  Then 


[x  < c] 


If  ^ 0 

“■"•vT  - 


/ \ C 

The  probability  density  fanction  o_(*)  of  X - F is  obtained  by 
differentiating  (2l): 


(21) 


-h2- 


q(\') 


1 r (x  <-  '}  XU;)’  y U'Y  ) “ ( ^ - '') 


(??) 


Wfion  the  ^npuL  to  the  cin{;l.c  ]dnc  procesooi-  ccnsis'v,s  of  noise  alone, 
X has  a K'ylci£yi  dlatrlLution: 


P(X)  = 


2 cxiJ  - 


2 C7- 


and  Bo: 


q(^)  = 


exp  ( ) d X 


2cr-‘ 


XO;) 


" ■ r^'2  ’ “ ‘ 


. (1  - .X-p  . ^ ) . _1  ( t . ch 

2 ex'  ^ -2  > 


A ^ 

rnc  characteristic  function  of  X*^  is  then; 

J (u)  = E exp  (i  u X ^)  = / exp  (i  u^)  q ( ^ ) d ^ 


1 t ^ 1 

1 - e:qo  { -x)  + ~A  exp 

2^^^  1 > 


2 c^' 


P (1  - 2iuc--^) 


vhcrc’ : 


C - ni.'ix 


(c,  . 


llio  r.vt;ra.C'3  nunbcr  of  fiilr.c  alarrr.s  ic  V ET.i:ation  (l^)  and  is 

ri'peat-  d ‘belcv: 


in 


n--i 


khJ-^ 

„ ) 


I’r  l^r(j)  > cj 

n 

1 - [r(j)>  c| 

__  ^ 

vhero  mv  P 


[l'(o)  >cj 


is  given  "by  Equation  (23) 


W’.en  the  target  is  present  X has  a Rice  distribution: 

f X ‘‘ 

p(x)  n - exp  ( - ”2  ' ) 1 (x) 

2 


end  for  this  case,  from  Equation  (22) j the  probability  density  function  of 

Ao  I . 

X = ^ is  given  by: 


q(^)  - 


— CXIJ  (• 


2 4 


2cr 


2~)  Io(x)  d X 


? cv 


2 fcxj. 


(- 


t 


["i  - a (.  n rJ-  )1  ?(^) 


cx: 


?Q  - 


(-  ) I (Jr)  u(S;-c^) 


-ir,- 


? 


A? 

t)i  ■ CiUii-acL-jrirtic  fo.’iot-Jon  of  X ir.  t>. -i 


(u)  r:  exp  (i  U ) q (^,  ) d ^ 


■f 


= 1 - Q , c/c>^)  + exp  - 


2 


o-~>  K 

• y iV(-fV  >_ 

k=0  ' r-i 


/ 2 / 

( 2 c'  ■ /c'  ) 


r-6  (k-r)]  (l  - 2iucv'^)^'^' 


fro.a  which  the  characteri r.tic  fanction  of  r*(j)  is  ohi-ainod  af 


= f(v) 


This  can  bo  inverted  directly  to  obtain  the  pro’o<ability  de'nnity  foncoion 

r(j): 

1 n 

i:; . . . >: 


a p 

PCS')  -X.  h’ r ■ 

n-0 


k =0 
n 


p „ ^k,-*  ...  +k„ 

(c^A)^  ^ 

k^/  ...  k^i  (2c^") 


-2(r,H...+  r_)  2, 

: ■*■  ^ exy>  - n /P 


n n =0 


Y - r,^  ...+  r„  + n-1 

(c:q.  ( r - n 0^)  ^ " , y>n  c-P 


n 


Tnori: 


id(j)  - ^ (.'l’  ' ■ <:/•-  ) • ■ • ]> 

r:  0 ..  h*  = 0 k ^'o 


k,  4 . . .^  k r,  -t . . .-I  r 

ml . y y " (i-c!Ay.l ! J. £V!) 

V I % I ^ ^ ~ --  \l  /_  vl 


VI  V I ^ 

i * ' • • n • ^ '-0 

n 


(Vr,)!  ... 


1 - I ( - 


C - n c 


r,  -I  . . .-t  r,  ,H  n 


^ f ^ . . . ^ + n - J ) 


vhere  C - r..'',x  ^C,  n j-  . Trie  prokab H ity  cr  detection,  is 

plotted  vs.  the  avorac;.-j  mr.r.bcr  of  false  alanns,  with  c as  a lara.ietor 

in  Fiaurc  6. 


I'grforr.ance  of  Trireshol d-Sqi:ere  Taw  J-rocescer  Und>r  UV7A  Conditions 
The  average  ni'vber  of  false  alan.is  under  these  conditions  is  {^Lvon  by 
Equation  (?i'0.  Tne  single  piri3  jiroccssor  output  vith  tarci-it  present  under 
ITir'TA  conditions  is  denot/jd  X and  has  pi'obability  density  function  (see 
Equation  (20)): 


Or>  CV^ 

= z z; 

m- 0 ti-  0 


-V  ?nH 

p (r.,n)  X exp  - 


I>et  X be  the  output  of  a tlires'noldi  113  d'  vice  v'nose  input  is  X.  Trie 

/^p  j, 

probability  density  function  of  X = is  obtained  fron  Equation  (2?)  and 
is  £iven  by: 


p(^  ) - IV  ^ ^ 

n t- 

^ i X-  V.  /5(r','0  t,  «i.  (--S-J  u{J,-/) 

ti-6  n-d  2 

A p 

Tiic  c>iaract(.-r.' cjLJc  functioa  of  X'  is  facn: 


(u)  r Y,  exp  (iu  X^)  = / o;n  (i\iV)  p(  t-  ) d 


OO  oo 


2/o 

P nil  c p er 

- , n 


ir.=  0 n 0 


m i 


c>'> 


E-0  n=0 


Vfe.4  ) 

^ 2 c:^ 


■ 4 rn--vT. 


n-r 


(-1U  t ) 


r+I 


The  charecieristic  function  of  r'(j)  vlll  then  be: 

. „ "iJ 

$p(u). 


^(u) 


which  can  be  inverted  uninj  I’air  431  Kefcrcnce  5 lo  yiej.d 


->:8- 


r(r).  / 0 

y.  6 


V-  -- 


>,:  s:  ...s::  >: 

ri  :0  t’  . 0 n,-0 

■1-  i.  K iC 


-])  ... 


n,  n, 


I.  -X 

T*  t ^ 1‘,  - O 

3 ‘ 


G(ii 


k’ 


(rj^+  ...+  Tj.  ^ k-l)J 


( Y ~ y c^) 


r^+  ...+  rj^  k-1 


t->q>  (-  2 )>  5 > k c 

2 o“ 


vnerc : 


c*o  c>c» 


,2 


f(i,,n)  I ( „)  _ 

ir.  -O  n=0  2 cT  n H 


0(n,r) 


Tn-r)i 


(X) 


n-r 


cjrp  - __ 


2 c- 


?(-n,n)  = (n<)  ^'(m,n) 

and  ^(n,n)  ic  defined  in  C-clion  3 3. 


-hn- 


'c 


^ 1 <f> 
} 0 


/ (c) 


c-rt 


i:>::  -4:5: 


rij^O  iK=0 


rv„0  ryO 


rj)  .../T(n^,  n.)  V---  r— 

f ^ ...2^  0(„^,ri)  .. 

'■k'O 


W.'  ’'k> 


? r,^  ...-i  r,  + k 

(2c^')^  ^ 


1 - I ( 


C).  - k c' 


r,.  +■ 


, r^f..,+rj.  -t  k - 1) 


vaerc  : 


C|^  = mx 


y.  c J . Tlie  probability  of  dotection,  P^(j),  ii 


plott-jd  vs.  the  avcracc'  na*cr  of  false  alams  in  PjGi.Lro  7 for  o"^  ’ = 1 


and  In  Fierro  8 for 


10 


~Yj~ 


iC).  n i t,y  (>r  ivjv.-cUon  iv.lor^^io:!  r>chr.;i.;o 

Vac  tecL  TiL'k  li.;;  [^(j)  3C  to  to  ccnputfj^l  fur  ti;o  {aL'a  cxLoaa Ior.3 

coiTc  rjv^r. ’ lac  to  tho  larcoat  valu:^r.  of  X toid  t>>_-  1.cr;t  ftjrict.ic'aG  for  other 

A 

p.''.lh  cxtar.r.icn'-:  rr-e  to  to  dvc'.i'od.  liie  nov  pi-ote-til  iX.y  oI  doiocti on  h^(j) 


bcconor. : 

yj) 


^d.-tocticn  afLcr  j {'hijr,  trii  • tore ’t  rath  in  not  du  .i;od^ 


IJov/: 


• I'y  ^ true  tra'cot  iat;i  Ir  rot  dii  .c<- 1 J. 

= Pjj(j)  • Pj,  -^truo  tort^ot  jath  is  not  du.nr-'cd^ 


{- 


P trur;  tercot  j ath  is  not  dunpod 


} 


~|  I P_  J true  tar^ot  path  is  not  dunp-^i  C'l'  P*  pi“f  }■ 

k=2  t J 

[l  - P,  norecfH-1  indopendont  R.oylo i c’r.  landcr.  vai  la't  Vo 


exceed  one  Rice  ro.ndom  variable  J.  J 


J-1 


H-1 


1 - 


^ (V)  (•'r  {•'=0  > Xl]  ) u - Pj  {Xo  > Xj] 

kye 


H-k-1 


J-1 


('^7} 


r 


-51- 


vlicro  I a'j  a Riylo  j dir.Lj  ibu'J  on  and  ban  a Rico  di  rJ,r Ibu'oi  on : 

X y ^ 

r.  (y  ) . - .. 

Pq'^'O'  ' P Ci.^i  - - ~ 

CJ"  2.r^- 


PiCX])  - 


Xl^  ^ 


fX 


q.  ( - 


2c- 


) I„(Xj) 


(28) 


I-C‘t  d'noLc  the  dirRritmtion  funcLlon  of  X^: 

X 

- IV  {Xo  <- 


= 1 - exp  - 


•X 


2 o- 


Conditioning  on  X,  ve  have: 


ro 

'[Xo>Xij  =J  -|^Xq>Xp  I X;j^  = Xj  Pp(X)  d X 

Oo 

- Po(X)j  p^(x)  dX 


2cr* 


exp 2 


I (X) 


cv 


2 cr- 


d X 


*=  -2'  exp  - — 


4“ 


.<■ 


Suhr. t. i iu‘.iri"  in  {?'()  vc  obt,:n 


•:i  1 n : 


Pj.  •/"true  tra'C'.-i  j;-.  not  dn.Ti;.f.‘d^ 


H-1 


’ - - -1-  )'  (1  - 


exp  - j 


V^_) 

Pr/’jT 


Under  UVT<\  condl tionr.,  X in  roplfieed  by  X end  Pj('^)  oj"  Fqunti  )i 
is  reiOaecd  by  (r.oc  Dquation  (l^l-))  : 


Ei-O  n-0 


2ntl  2 

n)  X exp  - 

2c:r^ 


L'ncn; 


^ rv 


Co 


cyO 


(exp  - — p ) ( ^ ^ exp  - dX 


2 

2 


CO 


/ \^  / 2vH!+n  / X 

(-1)  ( ex'  ) (l;  2;  rr>n) 


X I 

n-5  n=6  raj  n{  (nun)] 


= P. 


(29) 


U- 


if 


M 


SubstltuLinc  (29)  in  (?.'()  wc  obtain: 

Pj,  ^truo  target  is  not  dumped  under  UVd'A  conditions^ 
p H-1 

1 - 2-  J k n,.  ( 1 - Pf ) 

^D^J) 

VJ) 


'W- 


77ie  d.-tradition  factor  P^( J ) /Pjj( .! ) in  plotted  vn.y/./i!  v.'th 
an  a para..,:;tcr  for  additive  noire  corditionr  in  Fif^ur'^  9 arid  for 
UViYi.  conditionn  in  Fi£7.u\;-  10. 
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